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1.  INTRODUCTION 


The  propagation  of  neutral  particle  beams  in  the  upper  atmosphere  is  of 
substantial  scientific  interest.  Since  the  time  of  Alfven  [1 J  this  interest 
has  been  motivated  by  the  consideration  of  natural  phenomena  such  as  the  origin 
of  solar  systems  and  aurorae.  Most  recently  there  has  been  considerable 
interest  in  this  problem  because  of  the  possibility  of  artificially  created 
particle  beams. 

The  most  important  parameter  of  a  neutral  beam  is  determined  by  its 
energy  density.  Thus  it  is  imperative  to  understand  the  energy  loss  mechanisms 
of  a  neutral  beam  as  it  propagates  in  the  upper  atmosphere.  The  energy  loss 
mechanisms  can  be  divided  into  three  main  categories:  1)  dispersive  effects 
which  reduce  the  energy  of  the  beam  by  lowering  its  density.  2)  particle  energy 
loss  mechanisms  which  reduce  the  energy  per  particle  in  the  beam  and  3) 
stripping  processes  which  lower  the  neutral  beam  energy  density  by  ionizing 
particles  in  the  beam  and  thus  subjecting  them  to  the  influence  of  the 
geomagnetic  field. 

Depending  on  the  energy  and  current  density  of  the  neutral  beam,  there 
are  six  main  physical  processes  which  will  diminish  the  beam  energy  density: 

1)  the  geometric  dispersion  of  the  beam, 

2)  scattering  processes  due  to  weak  collisions  with  atmospheric 
constituents , 

3)  nuclear  interactions, 

4)  radiative  processes, 

5)  ionization  processes  due  to  strong  collisions  with  atmospheric 
constituents , 

6)  self-ionization  of  the  beam. 

The  first  process  is  an  effect  determined  by  the  generation  of  the  beam 


and  depends  on  the  degree  of  collimation  of  the  beam.  It  causes  no  particle 
energy  loss  and  only  reduces  the  density  of  the  beam  by  broadening  it.  The 
second  process,  scattering  collisions,  reduces  both  the  energy  per  particle  and 
density  of  the  beam  but  is  relatively  unimportant  for  high  energy  beam 
particles  since  they  would  become  ionized  before  undergoing  significant  energy 
loss  due  to  scattering  collisions  [2J.  The  third  and  fourth  processes  may  be 
neglected  for  neutral  beam  energies  less  than  30  MeV  [2J.  Obviously  these 
processes  would  become  much  more  significant  and  even  dominate  the  energy  loss 
mechanisms  for  high  enough  beam  energies.  The  fifth  and  sixth  processes 
provide  the  main  particle  energy  loss  mechanisms  in  the  1  MeV  energy  range  and 
therefore  are  the  most  effective  in  reducing  the  energy  density  of  the  neutral 
beam.  This  Report  examines  the  last  process  in  order  to  determine  its 
relative  importance  with  respect  to  other  energy  loss  mechanisms. 

2.  CROSS  SECTION  DATA  REVIEW 

When  a  neutral  hydrogen  beam  propagates  through  the  upper  atmosphere 
there  is  some  probability  for  collisions  with  the  atmospheric  constituents. 
Numerous  experiments  have  been  performed  over  the  years  to  investigate  the 
effects  of  atom-atom  and  other  collisions  among  the  atmospheric  constituents 
[3J.  These  experiments  have  determined  that  there  are  three  main  interactions: 
1)  stripping  or  separation  of  the  beam  hydrogen  atom  into  a  proton  and 
electron,  2)  ionization  of  the  atmospheric  or  target  particle  with  the 
resultant  release  of  a  free  electron  and  ton,  and  3)  electron-ion  (proton) 
recombination.  Most  of  the  data  available  in  the  literature  [3j  give  results 
for  the  total  cross  section  for  the  production  of  various  reaction  products 
over  a  wide  range  of  energies.  These  results  demonstrate  that  for  beam 
energies  above  LOO  KeV  the  cross  section  for  recombination  is  so  low  that 
this  process  may  be  neglected.  It  can  also  be  observed  that  the  cross  section 
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1  for  stripping  or  ionization  is  essentially  proportional  to  the  number  of 

electrons  in  an  atom  and  further  that  there  is  negligible  difference  in  the 
stripping  cross  section  due  to  a  neutral  atom  or  due  to  its  ion  and  electron 
components  [4J. 

The  more  recent  experimental  cross  section  investigations  are  worth 
noting  in  more  detail.  Current  investigations  are  being  conducted  mainly  by 
four  groups.  Gilbody  et  al  [5]  have  performed  numerous  experiments  and 
obtained  total  cross  section  results  over  a  wide  range  of  energies  for  proton 
beams  incident  on  various  gases.  This  group  has  measured  cross  sections  both 
for  ionization  and  excitation  of  gases  due  to  proton  bombardment.  Van  Zyl  et 
al  [bj  have  investigated  the  interaction  of  neutral  Hydrogen  with  Nitrogen  and 
Oxygen  molecules  but  their  reported  results  have  been  restricted  to  very  low 
beam  energies  for  application  to  auroral  phenomena.  Toburen  et  al  17  J  have 
investigated  proton-atom  collisions  over  a  wide  proton  beam  energy  range  and 
reported  differential  cross  section  results  for  molecular  Hydrogen  and 
Nitrogen.  Rudd  et  al  [8]  have  also  performed  a  number  of  proton-atom 
experiment^.  In  addition  they  have  reported  cross  section  results  differential 
in  energy  and  angle  for  the  electrons  emitted  by  stripping  or  ionization  of 
neutral  Hydrogen  in  Hydrogen-Helium  collisions  at  intermediate  energies.  These 
differential  cross  section  results  clearly  indicate  the  presence  of  two 
distinct  groups  of  electrons.  One  group  appears  with  very  low  velocities  and 
is  presumably  the  result  of  ionization  of  the  stationary  Helium  target.  The 
other  group  of  electrons  has  a  clearly  defined  velocity  distribution  centered 
around  the  equivalent  atomic  Hydrogen  beam  velocity  and  is  presumed  to  arise 
from  the  stripping  of  the  neutral  Hydrogen  atoms  in  the  beam. 

Unfortunately  no  detailed  differential  cross  section  results  are 
available  for  neutral  atomic  Hydrogen-Nitrogen  or  Hydrogen-Oxygen  collisions  in 


3 


the  1  MeV  or  greater  energy  range.  Also,  while  the  investigations  of  Rudd  et 
al  18 J  have  yielded  useful  information  about  the  electron  products  very  little 
is  known  experimentally  about  the  resulting  ion  distributions. 

3.  CHARGED  PARTICLE  DRIFT 

The  charged  particles  produced  from  the  neutral  hydrogen  beam  stripping 
and  atmospheric  ionization  processes  will  have  different  behaviors.  Most  of 
the  ions  and  electrons  produced  by  ionization  of  the  neutral  atmospheric  atoms 
can  be  assumed  to  gain  little  kinetic  energy  and  thus  will  essentially  remain 
stationary  and  contribute  to  the  ambient  plasma  background.  Since  the  combined 
electron-ion  effective  cross  section  for  stripping  of  the  neutral  beam  is 
practically  equivalent  to  that  of  the  neutral  atom  [ 4  J  their  effect  on  beam 
propagation  is  negligible.  The  fast  protons  and  i  lectrons  which  arise  from  the 
stripping  of  the  neutral  hydrogen  beam  by  the  atmospheric  constituents  cannot 
be  dismissed  in  a  similiar  manner.  Obviously  this  stripping  reduces  the 
density  of  the  neutral  beam.  Also,  numerous  recent  experiments  [9j  have 
indicated  that  a  sufficient  density  of  neutralized  ions  will  establish  a 
polarization  electric  field  and  propagate  transversely  to  a  magnetic  field  with 
a  transverse  velocity  almost  equal  to  the  initial  velocity  of  the  neutralized 
ions  (.protons). 

Charged  particle  motion  across  magnetic  field  lines  was  first  proposed 
by  George  Schmidt  [lOj.  His  theory  has  since  been  refined  by  a  number  of 
authors  (11 j  and  confirmed  by  various  experiments  [9],  Peter  and  Rostoker  [11 J 
have  established  a  criterion  for  the  initiation  of  charged  particle  drift 
accross  magnetic  field  lines.  The  condition  is  that  the  initial  neutralized 
ion  beam  energy  should  exceed  the  energy  density  necessary  to  establish  the 
driving  electric  field.  In  other  words  the  static  dielectric  constant  must  be 
greater  uian  the  neutralized  plasma  mass  ratio.  For  a  hydrogen  (proton)  beam 


we  have 


They  quote  experimental  evidence  to  demonstrate  that  the  polarization  drift 
will  occur  for  the  static  dielectric  constant  e  >  200.  This  condition  is 
satisfied  for  ion  (proton)  densities  n^  >  10^  part/cnP. 

The  theory  of  Peter  and  Rostoker  l 1 1 J  is  applicable  for  large 
gyroradius  neutralized  ion  beams  injected  from  a  field  free  region  to  a  region 
of  high  magnetic  field.  Large  gyroradius  beam  implies  that  the  ion  (proton) 
gyroradius  is  much  larger  than  the  neutralized  ion  beam  radius.  It  does  not 
appear  that  a  large  magnetic  field  is  a  necessary  component  of  the  theory.  A 
more  serious  defect  of  the  theory  is  the  tacit  assumption  of  the  formation  of 
the  polarization  charge  layers  which  initiate  the  plasma  drift.  Experimental 
evidence  has  conclusively  demonstrated  that  these  charge  layers  are  formed 
when  a  neutralized  beam  is  injected  into  a  vacuum.  However,  the  situation  is 
different  in  the  upper  atmosphere  (a)  because  of  the  presence  of  the  ambient 
plasma  background,  and  ( b)  because  of  the  large  ion  gy roradius-beam  radius 
ratio  prevailing  in  the  geomagnetic  field.  It  is  conceivable  that  the  ambient 
plasma  may  act  as  a  short  to  prevent  the  formation  of  a  virtual  anode  and  thus 
the  polarization  electric  field.  The  theory  of  Peter  and  Rostoker  [11]  should 
be  carefully  examined  to  investigate  this  possible  effect. 


4.  THE  EFFECT  OF  BEAM  INDUCED  STRIPPING 

From  the  experimental  cross  section  results  [8]  we  know  that  the 
velocity  of  the  electrons  produced  by  the  atmospheric  stripping  of  the  neutral 


beam  suffer  little  degradation  of  their  equivalent  initial  atom  velocity.  We 
assume  that  this  is  also  true  for  the  protons  produced  by  the  same  stripping 
reactions.  Once  a  critical  ion  (proton)  density  is  reached  the  stripped  ions 
and  electrons  will  propagate  across  the  magnetic  field  lines  with  the  neutral 
beam  at  some  small  relative  velocity  difference.  These  ions  and  electrons 
will  then  enhance  the  stripping  of  the  neutral  beam  and  therefore  diminish  its 
effective  propagation  distance.  The  importance  of  this  effect  can  be 
quantified  by  adopting  a  simple  model. 

In  order  to  investigate  the  relative  importance  of  the  ionization  of 
beam  particles  due  to  collisions  with  atmospheric  constituents  versus  self 
ionization  of  the  neutral  beam  we  have  made  a  number  of  simplifying 
assumptions.  First  we  ignore  beam  dispersion  and  scattering  processes  which 
lower  the  beam  density  by  broadening  the  beam.  Second  we  ignore  atmospheric 
density  variation  with  height  (this  would  be  strictly  applicable  to  horizontal 
beam  propagation).  Neither  of  these  assumptions  are  critical  and  can  be 
accounted  for  in  this  model  fairly  easily  but  both  would  increase  the 

complexity  of  the  numerical  calculations.  We  also  have  made  no  attempt  to 

provide  a  detailed  mechanism  which  allows  the  stripped  ions  and  electrons  to 

propagate  with  the  neutral  beam.  Clearly  for  the  case  of  neutral  beam 

propagation  parallel  to  the  magnetic  field  this  is  not  a  problem.  The 
previous  section  has  reviewed  relevant  references  and  experiments  which 
provide  a  mechanism  for  neutralized  ion  drift  transverse  to  a  magnetic  field 
but  no  detailed  theoretical  or  experimental  treatments  exist  for  ion 
propagation  at  an  arbitrary  angle  to  the  magnetic  field.  Our  simple  model 
assumes  that  this  will  occur  and  introduces  a  geometric  factor,  y>  as  an 
adjustable  parameter.  This  parameter,  y>  can  be  thought  of  as  the  ratio  of 
electrons  (ions)  travelling  within  the  beam  to  the  total  number  of  electrons 


(ions)  produced  by  the  stripping  process  and  is  allowed  to  vary  between  0  (no 
electrons  (ions)  travelling  with  the  beam)  and  1  (all  electrons  (ions) 
produced  by  the  stripping  process  confined  within  and  travelling  with  the 
neutral  beam). 

Under  these  assumptions  the  change  in  beam  density  with  distance  of 
propagation  is  given  by 

=  ~o[aa  +  a+a-f.  +  a_a_  +  y+0+b+  +  y_0_b_Jb  (4.1) 

where  x  is  the  propagation  distance 

a  is  the  stripping  cross  section 

a+  are  atmospheric  neutral,  ion  and  electron  densities 

oc+  are  cross  section  corrections  to  account  for  species  variations 

b+  are  beam  neutral,  ion  and  electron  densities 

tU  are  cross  section  corrections  to  account  for  species  variations 
y+  are  geometric  factor  corrections 
It  is  obvious  that  the  rate  of  production  of  beam  ions  and  electrons 
are  equal  and  equal  to  the  rate  of  loss  of  beam  neutrals.  Thus 


3h+  =  3b_  =  -3_b 
3  x  3  x  3  x 


(4.2) 


We  can  assume  that  the  beam  ions  and  electrons  are  confined  within  the  beam 
and  allow  the  geometric  factors  (y+)  to  account  for  the  density  variation. 
Then  (4.2)  becomes 


dhf  _  db_  _  -db 
dx  dx  dx 

which  can  be  integrated  immediately  to  yield 
b+  =  b_  =  N-b 

where  N  is  the  initial  density  of  the  beam. 


(4.3) 


(4.4) 


Clearly  (4.4)  is  not  true,  for  two  reasons.  First,  when  the  relative 


velocity  difference  between  the  ions  and  the  beam  is  not  ignored  then  b+  are  not 
equal  but  probably  less  than  N-b.  However,  this  is  the  limiting  (worst)  case 
and  the  discrepancy  can  be  accounted  for  by  the  geometric  factors,  y±*  Second, 
if  the  velocity  difference  between  the  neutral  and  ion  beams  is  not  Ignored, 
then  b+  are  not  limited  to  N-b  and  could  in  fact  exceed  this  value  at  some 
distance.  In  this  model  there  is  no  difference  between  distance  and  time. 

These  two  variables  are  equivalent  and  simply  related  by  the  beam  velocity. 

The  removal  of  this  equivalence  can  be  accomplished  by  replacing  equation 
(4.3)  above  by  a  true  divergence  equation  which  is  not  symmetric  in  distance 
and  time.  Such  a  refinement  is  presently  being  investigated  and  some  of  the 
aspects  are  discussed  below. 

Using  (4.4)  in  (4.1)  yields 

^  =  -olaa  +  a-|.a+  +  a_a_  +  (y+S+  +  Y-3-)  (N-b)  Jb  (4.5) 

which  can  be  written  as 

^  =  -(A+BN)  b  +  Bb2  (4.6) 

dx 

whe  re 
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b 

A-t-B(N-b) 


=  exp 


~(A+BN) (x+C) 


(4.9) 


With  the  initial  condition  that  at  x=o,  b=N  we  can  eliminate  C  from  (4.9). 
Simplifying  yields 


b  N 

A+B(N-b)  ‘  A  eXp 


-(A+BN)x 


(4.10) 


This  equation  (4.10)  specifies  the  beam  density  as  a  function  of  propagation 
distance  and  known  or  assumed  constants.  In  the  limit  of  no  beam-beam 
interaction  B+0  and  (4.10)  reduces  to 

b=N  exp  (-Ax)  (4.11) 

which  is  the  correct  behavior. 

From  (4.10)  it  is  easy  to  see  that  there  is  no  combination  of  constants 
which  leads  to  zero  beam  density  for  any  finite  propagation  distance.  It  is 
also  apparent  that  the  neutral  beam  density  decreases  more  rapidly  when  beam 
effects  are  included.  A  comparison  of  equations  (4.10)  and  (4.11)  leads  to 
the  conclusion  that  beam  self  ionization  effects  become  important  when  BN>A. 

If  we  ignore  minor  variations  in  the  cross  sections  and  geometric  factors  this 
implies  that  beam  effects  will  be  Important  when  the  density  of  the  beam  is 
equal  to  or  greater  than  the  atmospheric  density. 

This  can  be  Illustrated  by  substituting  the  appropriate  data  into 
equation  (4.10).  At  200  km  altitude  the  neutral  atmospheric  density  a  ~  lO^O 
part/cm^  [4J  while  the  maximum  ambient  plasma  densities  are  a+  ~  10^  part/cm^. 
Since  the  ambient  plasma  densities  are  at  least  four  orders  of  magnitude  lower 
than  the  neutral  density  they  can  be  safely  ignored.  At  a  beam  energy  of  1  MeV 
there  is  very  little  variation  in  the  cross  sections  for  the  various  processes. 
Thus  we  can  take  o  =  2x10“^  cm^/part  [  3  j  and  set  all  the  cross  section 
corrections  equal  to  1. 

We  have  chosen  to  let  the  beam  density  be  variable  and  given  by 
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N  ■  .446x10^  part/ctP  (times  beam  current  in  mA). 


These  approximations  have  been  made  only  to  simplify  the  calculations 
and  are  by  no  means  an  integral  part  of  the  derivation.  It  would  be  a  simple 
matter  to  use  more  precise  numbers  but  equation  (4.10)  is  relatively 
insensitive  to  small  changes  in  these  constants  and  more  precise  numbers  would 
not  significantly  alter  the  results. 

The  solution  of  equation  (4.10)  can  be  obtained  for  any  desired  beam 
density,  distance  and  altitude.  We  have  chosen  to  present  the  solution  in  the 
form  of  Figures  1,  2,  and  3.  These  figures  are  illustrative  of  the  strengths 
of  the  various  processes  and  are  useful  for  discerning  trends.  However,  in 
light  of  the  various  approximations  inserted  for  the  numerical  constants,  one 
should  understand  the  limitations  of  the  numerical  results.  In  other  words, 
these  figures  accurately  portray  the  shapes  and  relationships  of  the  curves 
(and  underlying  processes)  but  do  not  yield  accurate  numerical  results.  Figure 
l  shows  the  distance  from  launch  of  the  neutral  beam  where  a  given  beam  density 
ratio  of  .1  (survival  probability)  would  be  obtained  as  a  function  of  the 
geometric  factor  y  for  various  values  of  the  beam  current  (in  mA).  There  is  no 
beam  induced  stripping  when  y=0.  The  absence  of  beam  induced  stripping  would 
be  indicated  by  a  horizontal  line.  As  is  clear  from  the  Figure  for  neutral 
beam  currents  of  10  mA  or  less  beam  induced  stripping  effects  are  negligible  at 
200  km  altitude.  A  neutral  beam  current  of  lO^raA  yields  a  beam  density  which 
is  of  the  same  order  as  the  neutral  atmospheric  density  at  200  k.m  altitude.  It 
is  obvious  from  the  Figure  that  beam  induced  stripping  is  important  in  that 
case.  Figures  2  and  3  are  the  same  as  Figure  1  except  for  different  beam 
density  ratios  or  survivai  probabilities. 

It  should  be  emphasized  that  all  the  Figures  are  for  a  fixed  altitude. 
What  is  important  to  note  is  that  beam  induced  stripping  becomes  important  when 


10 


the  neutral  beam  density  is  of  the  same  order  as  the  atmospheric  density.  In 
other  words  for  fixed  altitude  beam  induced  stripping  increases  as  you  increase 
the  neutral  beam  density.  Also,  for  a  fixed  beam  density,  beam  induced 


4 


stripping  becomes  the  dominant  process  as  you  increase  the  launch  altitude. 

Finally,  we  wish  to  sketch  how  a  more  accurate  description  of  the 
electron  production,  with  relative  electron  motion  can  be  formulated. 

The  electrons  emerge  with  somewhat  less  energy  from  the  stripping 
process  than  what  they  had  as  constituents  of  the  neutral  beam  in 
the  beam  frame.  This  leads  to  an  average  backward  relative  velocity, 
say  w.  If  the  electron  source  function  at  point  x'  measured  from  the  front 
edge  of  the  beam  and  at  time  ^'is  S(x',t'),  the  electron  density  at  _x  and  £  is 
is  given  by 


b_(x, t)  “  /  S(x'  =,x-w[ t-t '  ] ,  t')  dt' 

t-x/w 


(4.12) 


The  lower  limit  tm  of  the  integration  is  fixed  by  the  condition  that  the  most 


distant  point  from  which  the  electrons  can  propagate  is  the  beam  front, 
x'  *  0,  i.e.  x-w  [t-t'm]  =  0 
On  the  other  hand,  obviously 


S(x, t) 


(x,  t) 


(4.13) 


(4.14) 


An  equation  for  the  time  evolution  of  b(x,t),  analogous  to  (4.1),  can  now  be 


written  down  (ion  effects  neglected  here) 


n  L 

—  =  -  o(w)  wlaa  +  a+a+  +  cx_a_  -hy-B-b-jb 

O  t 


(4.13) 


Combined  with  (4.12)  and  (4.14)  and  yields 


j-jr  (x,t)  =  -[A(w)-B(w)  /£_x/w  (x-w  [t-t'  J  ,t'  )  dt  j  b(x,  t)  (4.16) 


iff 


4  “  o(w)w  [aa  +  a+a+  +  a_a_J 
B  =  a (w)w  y_3_ 


(4.17) 


(4.16)  is  a  complicated  functional  differential  equation  whose  solution  is 
under  investigation. 


5.  GEOMETRIC  FACTOR  CALCULATION 

One  can  note  from  the  results  of  the  previous  section  that  a  lmA/cm^ 
beam  propagating  horizontally  at  200km  altitude  suffers  negligible  beam 
induced  stripping  compared  to  the  stripping  caused  by  the  atmosphere.  It  is 
clear  that  beam  induced  stripping  becomes  a  more  important  energy  loss 
mechanism  for  higher  beam  densities.  It  is  also  apparent  that  beam  Induced 
stripping  will  become  the  dominant  energy  loss  mechanism  at  higher  launch 
altitudes  since  the  neutral  atmospheric  density  falls  off  rapidly  with 
altitude.  Thus  we  have  attempted  to  gain  a  better  understanding  of  the 
geometric  factor  for  the  case  of  beam  propagation  parallel  to  the  magnetic 
field  by  adopting  a  simple  model  in  which  the  ions  stripped  from  the  beam 
continue  with  the  neutral  beam  and  the  stripped  electrons  form  a  neutralizing 
cloud  around  the  beam. 

In  cylindrical  coordinates  the  equation  of  motion  for  an  electron  is 


••  e  eB  •  *2 

r  =  —  e  +  —  rd>  +  r«t 
in  m 


2m  mr 


(5.1) 


where  E  is  the  self  induced  electric  field,  B  is  the  magnetic  field  (assumed 
constant),  r  and  <j>  are  polar  coordinates,  z  is  the  direction  of  the  beam  and  C 
is  a  constant  of  motion.  Energy  conservation  yields  the  first  integral  and 
after  simplification  (5.1)  can  be  written  as 


m 


Y  v 


(5.2) 


•2,1  (U-eV-wHC)-wHr2- 

m2r2 

where  U  is  the  total  energy,  U  and  C  are  constants,  V  is  the  scalar  potential 

,  eB 

and  «H  “  ^ 

From  the  diffusion  equation 


$.j  “  0 


(5.3) 


and  assuming  isotropy  in  and  z,  we  get 

7  t  (prr)  =  0  (5.4) 

r  or 

or 

prr  =  constant  =  Y  (5.5) 

Thus  the  electron  density  p  is  given  by 


Y 

P  *  T 
rr 

From  Poisson's  equation 
V2V  =  -4irp 
We  can  write 


3 V  -4ttY  ,  dr 
9r  =  r  J  r 


Therefore,  substituting  r  from  (5.2), 


_3V 

a 


V  -4ttY  f 

7--T  J  I 


dr 


C2  31/2 


T  2  2  &  1 

^[(U-eV-u,HC)-u)Hr  - 


(5.6) 


(5.7) 

(5.8) 


(5.9) 


which  is  an  equation  for  the  scalar  potential  solely  in  terras  of  r  (the  radial 
distance  from  the  center  of  the  beam)  and  constants.  Unfortunately  this 
equation  has  no  analytic  solution  and  ultimately,  it  must  be  solved 
numerically.  However,  an  insight  may  be  obtained  by  using  an  iterative 
technique.  We  first  assume  there  Is  negligible  external  potential  V(r).  We 
then  can  calculate  the  first-order  generated  potential  and  deviation  of  the 
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orbits. 


Using  this  approximation,  eq.  (5.9)  becomes 


3V  4irY 


lfv — 


HC)-^Hr;-  ~  2'J 


(5.10) 


which  can  be  integrated  to 


3V  4uY  .  _i  H 


3r  =  ojjjr 


-mw2r+U-wnC 


[U(U-2whC)J 


(5.11) 


where  Ej  is  an  integration  constant.  The  boundary  conditions  here  are  set  at 


the  turning  points  of  field-free  particle  trajectories: 


j  U-0)HC-  Ju(u-2uhc)  1 1 2 


U-u)HC-ty|U(U-2o)HC)j1/2 


(5.12a) 


(5.12b) 


Equation  (5.11)  gives  the  electric  field  due  to  the  electrons  only.  We 
expect  this  to  vanish  at  the  lower  limit  rj .  This  boundary  condition  yields 

-rau)2r  +  U-uqjC  -i 

3 V  _  4n Y  .  H _  .  2tt2Y  ... 

3  r  ~  wHr  sln  [  [U(U-2whC)J  1/2 J  uHr!  (5.13) 

The  constant  Y  may  be  evaluated  from  consideration  of  the  total 

electric  field,  due  to  the  electrons  and  ions.  The  field  due  to  the  ions  may 


be  written  as 


2-nrp^  inside  the  beam  (r<rQ) 
2ur^p  outside  the  beam  (r>rQ) 


(5.14) 


where  we  assume  a  uniform  ion  charge  density  pi  inside  the  beam  which  has 


radius  r0. 


The  total  electric  field  is  just  the  sum  of  (5.13)  and  (5.14).  Using 


Che  boundary  condition  that  the  total  electric  field  is  zero  outside  the 
charge  region  (r  >  r2) ,  we  get 


2 

,  ^Hr°pl 

w(l  +  r2/  ri ) 

We  can  now  write  the  expression  for  the  total  electric  field: 


(5.15) 


_  4roP  1  .  ,  2  2.  .  ,  .  rr  )  2ir r0p i 

Etotal  r  l-i-  r2/r1  sin  [U( 0-2^0]!/ 2  (r2_r  ^  1  2  J  rt  +  r2 


1  [_  J  I  (5.16) 

Unfortunately,  (5.16)  cannot  be  analytically  evaluated  to  give  the 
potential  V(r).  Thus,  we  approximate  the  potential  which  generates  (5.16)  by 


a  parabolic  form 


Ve(r)  “  ar2  +  br  +  c  for  r>rQ 


VA(r)  =  -ar2  +  6 
with  boundary  conditions: 


for  r<rr 


(5.16a) 

(5.16b) 


vi ( r0)  -  Ve(r0)  (5.17) 

v[(r0)  =  Vg(r0) 

Ve(r2)  -  0 

Vg(r2)  *  0 

V2Vi(r)  *  -4trp r<rQ 

These  boundary  conditions  allow  determination  of  the  5  coefficients  in 
equations  (5.16).  The  result  is 


V(r)  --^2-  (r-r22) 

r2  ro 

=  -a(r2-r0r2) 


(5.18) 


where 


a  -  irpj. 

We  now  use  this  potential  to  solve  for  the  new  turning  points  of  the 
motion.  This  yields  an  estimate  for  the  geometric  factor  y 
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It  is  illuminating  to  write  the  equations  for  rj  and  T2  in  terms  of 


characteristic  length  parameters, 


JL 

li  ea 


rL  = 


K 


2  _  c2 

^  mK 


(5.19) 


where  K  =  electron  kinetic  energy  at  the  point  of  scattering  (rQ)  and  the 
total  energy  U  is  given  by  U  =  K  +  ar(r0-r2).  In  terms  of  these  parameters 
and  equation  (5.18)  for  the  potential,  the  equations  for  the  turning  points  r2 
and  r,  are: 


2  2 

Rc  {  r0(ro-r2)  r2  Rc 

Rl  Re2  2rL2  2r2 


(5.20a) 


f  2  2,  ,  d2  Ks  .K  ,  2  n  n  (5.20b 

2  (r2-r,)  +  Rc  (2  -  2)  +  2  (ri-rQr)  =  0 

rL  4  ^  r2  rl  rE  1  ° 

These  equations  have  been  solved  numerically,  using  a  perturbation 
technique  with  initial  values  from  the  field-free  turning  points.  While  the 
perturbed  value  of  r2  can  vary  strongly  at  low  energies,  rj  is  not  very 
sensitive  to  these  parameters. 


6.  SUMMARY 

Present  research  efforts  can  be  divided  into  four  main  areas: 
l)  Investigation  of  "worst  case"  theory  of  beam  induced  stripping  to  obtain 
results  as  a  function  of  altitude  and  beam  energy.  2)  Refinement  of  "worst 
case"  theory  of  beam  induced  stripping  effects  by  obtaining  a  numerical 
solution  of  the  equations  which  break  the  distance-time  symmetry  of  the 
present  theory  and  explicity  account  for  the  velocity  variation  between  the 


beam  ions  and  the  neutral  beam.  3)  Calculation  of  the  geometric  factor  by 
solving  the  equations  for  the  electron  density  self  consistently.  4)  Investi¬ 
gation  of  the  polarization  drift  theory  to  account  for  the  presence  of  the 
ambient  plasma  background  in  the  atmosphere  and  for  the  possibility  of  the 
buildup  of  sufficient  neutralized  ion  density  to  initiate  the  polarization 
drift. 

7.  CONCLUSIONS 

We  have  demonstrated  through  a  simple  model  that  there  is  no 
catastrophic  self-induced  quenching  of  the  neutral  beam.  It  is  likely  that 
this  feature  will  remain  true  even  when  the  model  is  modified  to  include  time 
dependent  effects. 

We  have  shown  that  for  a  1mA/ cm^  beam  current  and  200  km  atmospheric 
altitude  beam  induced  stripping  is  an  insignificant  energy  loss  mechanism  for 
the  neutral  beam  compared  to  the  stripping  induced  by  the  atmospheric 
constituents . 

On  the  other  hand,  we  have  also  demonstrated  that  beam  induced 
stripping  is  significant  when  the  density  of  the  neutral  beam  and  the 
atmosphere  at  the  launch  altitude  are  comparable.  This  would  be  the  case  for 
higher  beam  densities  (currents)  at  a  given  altitude  or  for  a  given  density 
beam  at  higher  altitudes.  It  is  also  clear  for  high  launch  altitudes  that 
beam  induced  stripping  would  be  the  dominant  energy  loss  mechanism  for  the 
neutral  beam  since  the  neutral  atmospheric  density  decreases  rapidly  with 
increasing  altitude.  In  absolute  terms,  at  high  beam  current  densities 
(~  lA/cm^)  the  beam  induced  stripping  might  set  severe  limitations  on  the 
useful  propagation  range  of  the  beam. 
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8.  APPENDIX 


This  appendix  disucsses  the  application  of  the  theory  described  in  this 
Report  to  the  proposed  BERT  II  experiment.  It  is  clear  from  the  results 
obtained  in  this  study  that  the  stripping  (and  thus  the  energy  loss)  of  a 
IraA/ cm2  neutral  beam  due  to  beam  induced  effects  is  negligible  compared  to  the 
stripping  due  to  the  atmosphere  at  a  launch  altitude  of  200  km  and  thus  will 
not  be  an  important  process  in  the  proposed  BERT  II  experiment. 

The  stripping  of  the  neutral  hydrogen  beam  by  the  atmospheric 
constituents  and  by  the  stripped  beam  ions  and  electrons  both  cause  particle 
energy  loss  and  a  density  depletion  of  the  neutral  beam.  However,  the 
polarization  drift  may  cause  the  stripped  ions  and  electrons  from  both 
processes  Co  propagate  with  the  beam.  Thus  the  stripping  process  may  not 
diminish  the  total  energy  deposited  on  the  target  as  much  as  a  measurement  of 
the  neutral  beam  would  indicate.  It  would  be  very  useful  in  the  BERT  II 
experiment  to  be  able  to  distinguish  between  neutrals  and  ions  striking  the 
Target  and  to  measure  their  energy  and  velocity  distributions.  It  would  also 
be  important  to  measure  the  surface  potential  of  the  target  since  a  high 
charge  buildup  would  alter  the  ion  trajectories. 

As  the  particle  accelerator  being  developed  for  the  proposed  BERT  II 
experiment  is  not  100%  efficient  (with  a  relative  particle  yield  of  60% 
neutral  H,  20%  H+,  20%  H-  and  1 00%  e-)  the  theory  described  in  this  Report 
should  be  investigated  to  account  for  the  effect  of  the  ions  released  from  the 
accelerator.  One  should  also  carefully  examine  the  effect  of  the  increase  in 
neutral  atmospheric  density  around  the  accelerator  caused  by  rocket 
outgassing. 

Since  the  polarization  drift  effect  serves  to  enhance  the  total  energy 
delivered  to  the  target  it  is  important  that  it  be  thoroughly  investigated. 
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FIGURE  CAPTIONS 


FIGURE  1.  Plot  of  the  distance  from  neutral  beam  launch  (in  km  at  which  a 

given  beam  density  ratio  b/N  “  .1  (survival  probability)  would  be 
observed  as  a  function  of  the  geometric  factor  y  for  various  values 
of  the  beam  current  (in  mA/cm^). 

FIGURE  2.  Same  as  Figure  1  except  for  a  beam  density  ratio  b/N  =  .3. 

Figure  3.  Same  as  Figure  1  except  for  a  beam  density  ratio  b/N  =  .3. 
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